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Confirmation of persistent radio jets in the microquasar LS 5039 
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Abstract. We present here new observa tions conducted with the EVN and MERLIN of the persistent microquasar 
LS 5039 discovered by Paredes et al. (2000) with the VLBA. The new observations confirm the presence of an 
asymmetric two-sided jet reaching up to ~ 1000 AU on the longest jet arm. The results suggest a bending of 
the jets with increasing distance from the core and/or precession. The origin and location of the high-energy 
gamma-ray emission associated with the system is discussed and an estimate of the magnetic field at the base of 
the jet given. Our results suggest a well coUimated radio jet. We also comment on new observing strategies to be 
used with satellites and forthcoming detectors, since this persistent source appears to be a rather good laboratory 
to explore the accretion/ejection processes taking place near compact objects. 
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1. Introduction 

Microquasars are radio emitting X-ray binaries with rela- 
tivistic jets, like SS 433, GRS 1915-^105, GRO J1655-40 
or Cygnus X-3. The reader is referred to Mirabel & 
Rodriguez (199£) for a detailed review and to Castro- 



Tirado et al. ( |200lD and Durouchoux et al. ( p002[ ) for the 
most recent studies of microquasars. 

The V = 11.2 star LS 5039, located at an esti- 
mated distance of ^ 3 kpc and close to the galactic 
plane {I = 16.88 °, b — —1.29°), was initially proposed 



by Motch et al. ( |1997| ) to be the optical counterpart of 
the X-ray source RX J1826. 2-1450, a hkely High Mass 
X-ray Binary (HMXB). Using data from the NRAO VLA 
Sky S urvey and follow-up VLA observations, Marti et al. 
( |l998|) discovered that the source was also a non-thermal 
radio emitter with moderate variability. X-ray observa- 
tions of RX J1826. 2-1450 carried out by Ribo et al. ( [1999D 
showed that the X-ray spectrum was significantly hard up 
to 30 keV, with a strong Gaussian iron line at 6.6 keV, and 
neither pulsed nor periodic emission was found on time 
scales of 0.0 2-2000 s and 2-200 days, respectively. Paredes 
et al. ( ^000 ) discovered that the system displays relativis- 
tic radio jets at milliarcsecond (mas) scales, revealing the 
microquasar nature of LS 5039, and proposed an associa- 
tion with the 7-ray source 3EG J1824-1514. While in the 
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past the mass donor had been classified as an 07V((f)) 
star, optical and near-IR spectroscopic observations by 
Clark et al. ( ^OOlj ) sh ow th at it is in fact an 06.5V((f)) 
star. McSwain et al. ( ^001 ) obtained the radial velocity 
curve of the system, with their fitted parameters being 
a short period of P = 4.117 days, a high eccentricity of 
e = 0.41, a radial velocity of the system of Vb = 4.6 km s~^ 
and a mass function of f{m) — 0.00103 ^^0- Ribo et al. 
(2002) obtained optical and radio proper motions, and 
found that LS 5039 is a runaway X-ray binary with a space 
velocity ^ 150 km s~^. McSwain & Gies ( p002 ), based on 
wind accretion models, have recently suggested a neutron 
star as the compact object in LS 5039, and proposed an 
inclination of i ~ 30° for this binary system. 

The persistent radio emission of LS 5039 suggests that 
the jet is always present. With the aim to confirm this hy- 
pothesis, and also to detect the jet at larger angular scales 
than those imaged with the previous Very Long Baseline 
Array (VLBA) measurements, we observed this micro- 
quasar with the European VLBI Network (EVN) and 
the Multi-Element Radio-Linked Interferometer Network 
(MERLIN). In this Letter we present the obtained results. 



2. Observations and results 

We observed LS 5039 simultaneously with the EVN 
and MERLIN on March 1 2000 (3:20-7:10 UT, MJD 
51604.2) at 5 GHz. The EVN observations were performed 
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Fig. 1. EVN self-calibrated image of LS 5039 at 5 GHz 
obtained on March 1 2000. Axes units are in mas. The 
synthesized beam, plotted in the lower left corner, has a 
size of 7.60x6.96 mas in P.A. of —14°. The first contom' 
corresponds to 0.3 mJy beam^^, while consecutive ones 
scale with 3^^^. 
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Fig. 2. MERLIN self-calibrated image of LS 5039 at 
5 GHz obtained on March 1 2000. Axes units are in mas. 
The convolving circular beam, plotted in the lower left 
corner, has a diameter of 81 mas. The first contour corre- 
sponds to a fiux density of 1 mJy beam^^, while consecu- 
tive ones scale with 3^/^. 



Table 1. Flux densities at 5 GHz, length and P.A. of 
the LS 5039 structure features as detected by EVN and 
MERLIN. 
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MERLIN 




S 5 GHz 
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[mJy] 
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n 


[mJy] 


[mas] n 


Core 


29.3 






31.6 




NW jet 


2.6 


24 


-42 


4.0 


128 -29 


SE jet 


3.3 


34 


140 


4.2 


174 150 



Note: the errors in flux density, lenght and P.A. are 
0.1 mJy beam^^, 2 mas and 4° for the EVN features and 
0.4 mJy beam"\ 12 mas and 5° for the MERLIN ones. 



Inspection of the shortest MERLIN baselines reveals a 
constant flux density during the full run. That allowed us 
to image the source directly without splitting the data in 
time blocks. We present the obtained EVN and MERLIN 
images in Figs. and respectively. The MERLIN image 
is presented convolved with a circular beam, equivalent 
in solid angle to the interferometric synthesized beam of 
142x46 mas at position angle (P.A.) —47°. Our images 
clearly show that LS 5039 has a bipolar jet emanating 
from a central core. In both images the southeast (SE) 
jet is brighter and larger than the northwest (NW) one, 
as can be seen in Table |l|, where we have quoted the flux 
densities of the core and the jets and the lengths and P.A. 
of the jets. Their total lengths are ~ 60 mas in the EVN 
image and ^ 300 mas in the MERLIN one. 



with 7 antennas, namely (antenna, location, diameter) 
Effelsberg, Germany, 100 m; Jodrell Bank, UK, 25 m; 
Cambridge, UK, 32 m; Westerbork, The Netherlands, 
14x25 m; Medicina, Italy, 32 m; Nolo, Italy, 32 m; and 
Toruh, Poland, 32 m. Data were recorded in MklV mode 
with 2-bit sampling at 256 Mbps at left hand circular po- 
larization, yielding a full bandwidth of 64 MHz. The data 
were correlated at the Joint Institute for VLBI in Europe 
(JIVE). Standard fringe fitting and imaging analysis were 
carried out using AiPS and difmap. MERLIN recorded 
data with 2-bit sampling at dual polarization and a total 
32 MHz bandwidth. We analyzed the left hand circular 
polarization data excluding 1 MHz channel at both edges 
of the band, yielding a final bandwidth of 14 MHz. 



3. Jet parameters and intrinsic sizes 

These results confirm the existence of a two-sided radio jet 
in LS 5039 reported in pre vious VLBA 5 GHz observations 
by Paredes et al. ( 2000 ). This source does not present 
any strong outburst or, at least, none has been detected 
during the eleven month monitoring carried out by the 
Green Bank Interferometer between 1998 September 16 



and 1999 August 22 (Clark et al. gOOl|). On the other hand, 
inspection of the RXTE All Sky Monitor data (Levine 
et al. 1996) reveals that the X-ray flux was at the typical 
low level (Ribo et al. 1999| ) at the epoch of the EVN and 
MERLIN observations and several weeks before. All this 
suggests that the jets are persistent, as the VLBI images 
obtained up to now seem to indicate. 
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In the observations reported here, the jets extend fur- 
ther away than the ~ 6 mas of the VLBA observations at 
the same frequency. In all the images the jets have similar 
position angles, ^ 125° in the VLBA image from Paredes 
et al. ( |2000D , - 140° in our EVN image, and - 150° in 
our MERLIN image. These results suggest a bending of 
the jets with increasing distance from the core and/or pre- 
cession. 

The brightness and length asymmetry of the jet com- 
ponents may involve special relativity effects (Mirabel & 
Rodriguez 1999). Hence, assuming that this is the reason 
for the detected length asymmetry of the jets, we can es- 
timate some parameters by using the following equation: 



/3cos( 



(1) 



where p is the velocity of the jet flow in units of the speed 
of light, 9 is the angle between the direction of motion 
of the ejecta and the line of sight and /Xa and /Xp are 
the proper motions of the approaching (SE) and reced- 
ing (NW) components, respectively. Although we do not 
know the epoch of ejection of the terminal plasma, we can 
cancel the time variable by using the relative distances to 
the core c?a and dr, as expressed in Eq. |l]. We will iden- 
tify da and dr with the lengths of the approaching and 
receding jet, respectively. Using Eq. |l| and the EVN val- 
ues in Table ^ we obtain f3 cos = 0.17 ± 0.05, and hence 
13 > 0.17 ±0.05 and e < 80° ±3°. For the MERLIN values 
we obtain ^ > 0.15 ± 0.06 and 6* < 81° ± 3°. These values 
are similar to those previously derived from the VLBA 
image, of /? > 0.15 ± 0.04 and 6* < 81° ± 2°. We have not 
considered an eventual larger size for the NW jet with a 
brightness below the image noise level. McSwain & Gics 
(2002) have recently proposed an inclination oii ~ 30° for 
LS 5039. If we assume that the jet is perpendicular to the 
accretion disk, and that the disk lies in the orbital plane 
of the binary system, then 9 = i = 30°, and using the 
values from the EVN image, we obtain /3 = 0.20 ± 0.06, 
which indicates a mildly relativistic jet. 

The total length of the EVN and MERLIN jets is ~ 60 
and ^ 300 mas, respectively. Consi dering that the source 
is located at 2.9 kpc (Ribo et al. 2002), these angular 



lengths translate into linear lengths in the plane of the 
sky of ~ 175 and ~ 870 AU, respectively. Assuming that 
9 ~ 30° we obtain intrinsic total lengths of 350 and 
^1740 AU, respectively, and lengths of the SE jet arm 
of ^ 200 and ^1000 AU, respectively. Moreover, the jet 
width is smaller than one synthesized beam even in the 
EVN image. This implies a jet half opening angle < 6°. 

4. A scenario based on the 7-ray/radio emission 

It is interesting to use the presence of 5 GHz synchrotron 
emitting particles so far away in order to constrain the jet 
physical parameters. The starting point of the following 
calculations will be the scenario for the proposed EGRET 
emission of LS 5039 presented in Paredes et al. ( pOOOj ). 
In this context, the high energy 7-ray photons arise from 



inverse Compton (IC) scattering of UV photons from the 
luminous companion star by the relativistic electrons that 
later, after having lost part of their original energy due to 
IC losses, will account for the radio emission in the jets. 

We will adopt here a very simple model of an expand- 
ing jet. Cylindrical coordinates, z and r measured parallel 
and perpendicular to the jet axis, are the best choice in our 
case. The jet is assumed to form at a distance zq from the 
compact object and flows with a velocity v = (3 c. The lat- 
eral expansion of the jet is parameterized as r = rQ{z/zQY, 
where ro is the initial jet radius. A freely expanding con- 
ical jet would correspond to e = 1, while for slowed lat- 
eral expansion we would have e < 1 (see e.g. Hjellming 
& Johnston, 1988). The expansion velocity perpendicular 
to the jet axis is thus v±_ = dr/dt = tvr/z. Concerning 
the magnetic field, we will only consider the component 
perpendicular to the jet axis, because it has the slowest 
decay if conservation of the magnetic flux is assumed. The 
magnetic field along the jet will be thus parameterized as 
B = Bo{r/ro)-' = Bo(^/^o)-^ 

A relativistic electron injected at the base of the jet will 
decrease its energy mainly through adiabatic expansion, 
IC and synchrotron losses according to: 



dE _ 2 v^E 
'dt ~ ^3 r 



aicUr^dE' - asB'E 



(2) 



The factor 2/3 in the adiabatic expansion term comes 
from the lateral expansion of the jet. The constants aic = 
3.97 X 10^2 and as = 2.37 x 10~^ (c.g.s. units) are the co- 
efficients of the terms accounting for IC and synchrotron 
losses, respectively. The radiation energy density C/rad is 
assumed to be dominated by stellar UV photons from the 
06.5V star, whose luminosity above 10 eV amounts to 



~ 5 X 10^^ erg s ^ and J/j-ad = i*/[47rc(c 



2)]. With 



the recently determined orbital parameters (McSwain 
et al. |2001|, McSwain & Gies p003), the likely semima- 



jor axis of the orbit is a = 2.6 x 10 cm. Adopting these 
parameters, the radiation energy density close to the com- 
pact object is so high that electrons with energies higher 
than 10"'^ erg will easily produce 7-ray photons with ener- 
gies of 100-1000 MeV, as detected by EGRET. Moreover, 
integration of Eq. ^ taking only into account the IC losses, 
reveals that the energy of the electrons decays very fast to 
values of ~ 5 X lO"** erg, regardless of their initial energy. 

The IC contribution in Eq. || should be comparable to 
adiabatic losses at least at the base of the jet. Later on, 
the adiabatic expansion term will soon become dominant 
due to its slower decay (cx z^^). Therefore, at the base of 
the jet we can request: 



aic 



47rc(a2 + z§) 



2e vEq 
3 zq 



(3) 



We can use the values quoted above, together with v = 
0.2c and an energy of i?o = 5 x 10~^ erg for the elec- 
trons injected at the base of the jet, to solve for zq in this 
quadratic equation. The obtained results are 0.35, 0.85 
and 1.30 AU for e = 1, 1/2 and 1/3, respectively. The 
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synchrotron term can be shown to be negligible a posteri- 
ori and has not been considered here. 

Hereafter, we will proceed taking into account only the 
adiabatic expansion term, because it is dominant com- 
pared to the IC one, and rewrite Eq. || as: 



dE ^ 2e vE 
dt 3 z 



(4) 



This can be easily integrated to give E ~ Efj{z/ za)~'^'^^^ . 
As the energy decays, the synchrotron emission of rela- 
tivistic electrons will shift towards lower frequencies. In 
particular, the synchrotron frequency expected at a dis- 
tance z is given by: 



Vc = 6.27 X lO^^BE"^ 



6.27X 10i**Bo£^o(z/zo)-^^/^ . (5) 



The electrons still producing 6 cm (fc = 5 x 10^ Hz) syn- 
chrotron emission at about z = 1000 AU (the size of the 
MERLIN jets) are probably those originally injected with 
energies of i?o — 5 x lO^'* erg. Then, the magnetic field 
Bq at the base of the jet can be roughly estimated as: 



.0 x IQ-^^E, 



-1 7e/3 



(6) 



This provides values of 3.6 x 10^, 12 and 0.6 G for e = 1, 
1/2 and 1 /3, respectively. If magnetic flux is conserved, the 
corresponding field at the end of the MERLIN jet would 
be iJioooAU — 128, 0.4 and 0.06 G. For comparison, the 
only independent estimates of the magnetic field in the 
LS 5039 jets come from simple equipartition arguments. 
The results are merely indicative of the average magnetic 
field in the solid angle of the sky covered by the jets. Using 
the total fiux density and size of the MERLIN jets given 
m Table |l| we derive a field value of > 8 x 10 G, al- 
though this is only a lower limit because the jet w idth i s 
unresolved. The VLBA image by Paredes et al. (|2000|) , 
allowed them to estimate an approximate magnetic field 
of ~ 0.2 G for radio jets reaching up to a few mas. It 
is clear that a field of such intensity is more consistent 
with the e < 1 results than with the too high magnetic 
fields implied by a simple conical jet (e — 1). Therefore, a 
slowly expanding jet, and hence well collimated, is favored 
to power the LS 5039 radio emission up to ~ 1000 AU 
by the electrons that previously contributed to the 7-ray 
emission, without in situ acceleration being required. This 
is also in agreement with our upper limit of < 6° for the 
jet half opening angle. 



star exhibits intrinsic variations in the UV photon flux, a 
variability in 7-rays would be seen correlated to it. 

Since LS 5039 seems to have persistent radio jets, it 
would be interesting to study the possibility of detecting 
Doppler-shifted lines. As the jet flow has an estimated ve- 
locity of ^ 0.2c, the amount of Doppler-shift is expected 
to be comparable to that found in SS 433 (Kotani et al. 
[1996 ; Marshall et al. 2002). Hence, spectroscopic X-ray ob- 
servations with XMM or Chandra could be very useful to 
constrain the physics of the relativistic flow and to study 
the matter content within the jet. Finally, as pointed out 
by Levinson & Waxman (2001) and by Distefano et al. 



(2002), if the jets are hadronic we would expect the for- 
mation of TeV neutrinos, that could be detected in the 
future with detectors of km^-scale effective area. 

Overall, these results conflrm the existence of persis- 
tent relativistic radio jets in LS 5039, indicating that this 
source is a very good target to be studied with new in- 
struments in order to gain knowledge on jet physics. 
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5. Future prospects 

The high eccentricity of LS 5039 provides both, a variable 
accretion rate along the orbit implying a variable rate of 
electrons injected into the jet, and a variable radiation 
energy density close to the compact object. Hence, it is 
predicted a variability in the 7-ray luminosity correlated 
with the orbital period, which may hopefully be detected 
by INTEGRAL and/or GLAST. On the other hand, if the 
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